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Abstract 


*The  Korteweg- deVries  (KdV)  equation  and  the  finite-depth  equation  of 
Joseph  (1977)  and  Kubota,  Ko  and  Dobbs  (1978)  both  describe  the  evolution  of 
long  internal  waves  of  small  but  finite  amplitude,  propagating  in  one 
direction.'  In  this  paper,  both  theories  are  tested  experimentally  by 
comparing  measured  and  theoretical  soliton  shapes.  The  KdV  equation 
predicts  the  shapes  of  our  measured  solitons  with  remarkable  accuracy,  much 
better  than  does  the  finite-depth  equation.  When  carried  to  second-order, 
the  finite-depth  theory  becomes  about  as  accurate  as  (first-order)  KdV 

theory  for  oui\  experiments.  However,  second- order  corrections  to  the 

\  ' 

finite-depth  theory  also  identify  a  bound  on  the  range  of  validity  of  that 

) 

entire  expansion.  This  range  turns  out  to  be  rather  small;  it  includes 
only  about  hali  of  the  experiments  reported  by  Koop  and  Butler  (1981). 
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1.  Introduction 


y 


.^L. .. 


This  Is  the  first  in  a  series  of  three  papers  on  the  evolution  of  long 
internal  waves  in  a  stably  stratified  fluid.  If  these  waves  have  small 
amplitudes,  then  their  motion  is  governed  approximately  by  a  linear  wave 
equation,  with  small  but  cumulative  corrections  due  to  weak  nonlinearity, 
dispersion  and  dissipation,  and  possibly  to  a  slowly  varying  background. 
Several  theoretical  models  exist  which  Include  various  combinations  of  these 
cumulative  effects.  The  purpose  in  this  series  is  to  test  experimentally 
some  of  these  theoretical  models  of  the  evolution  of  long  internal  waves,  in 
order  to  obtain  some  notion  of  their  accuracy  and  range  of  validity. 

Our  experimental  program  was  similar  to  that  used  by  Hammack  and  Segur 
(1974)  to  test  the  KdV  equation  as  a  model  Of  the  evolution  of  long  surface 
waves.  In  fact,  these  experiments  were  performed  in  the  same  wave  tank. 
However,  in  this  series  we  are  less  concerned  with  a  comprehensive  test  of 
the  aspects  of  KdV  theory,  and  more  concerned  with  how  the  different  models 
fit  together  to  explain  the  variety  of  phenomena  observed  as  the 
stratification  is  varied. 

The  experiments  were  performed  In  a  wave  tank  which  contained  a  layer 
of  fresh  water  over  a  layer  of  brine,  and  bounded  above  by  a  free  surface. 
Some  representative  vertical  distributions  of  density  in  the  absence  of 
fluid  motion  are  shown  in  Figure  1  of  Part  III.  Because  the  density  varied 
continuously  with  depth,  the  linear  theory  predicts  an  infinite  set  of 
internal  wave  modes.  However,  we  are  primarily  interested  in  the  lowest 
internal  wave  mode,  so  we  will  consider  a  two  layer  model  with  a  free 
surface,  as  shown  in  Figure  1  of  this  paper.  This  simplified  model  admits 
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only  surface  waves  and  one  set  of  internal  waves. 

In  the  present  paper,  we  will  consider  two  models  which  are  weakly 
nonlinear  and  weakly  dispersive:  the  Korteweg-deVries  (KdV)  equation, 

fT  ♦  6ffx  ♦  fxxx  ■  0  (1) 

and  an  equation  due  to  Joseph  (1977)  and  to  Kubota,  Ko  and  Dobbs  (1978), 

♦  ffx  ♦  TCfxx]  *  0  •  <2,) 

where 

m 

T[f]  -  •f  coth|(x-y)  f(y)  dy  .  (2t>) 

and  the  integral  is  evaluated  In  the  principal-value  sense.  The  latter 
equation,  which  we  will  call  the  "finite-depth  equation",  may  be  written  in 
a  variety  of  equivalent  ways,  of  which  (2)  is  perhaps  the  simplest.  This 
list  of  equations  could  logically  Include  an  equation  proposed  by  Benjamin 
(1967)  and  later  derived  by  Ono  (1975) 

fT  *  ffx  -  HCfJ  •  0  .  (3) 

where  H[»]  Is  the  Hilbert  transform.  However,  we  will  not  consider  (3) 
because  we  have  no  experimental  data  in  the  range  of  parameters  where  (3)  is 


•  * ' 


< 


f 


val id. 

It  is  well-known  that  the  KdV  equation  describes  the  slow  evolution  of 
internal  waves  of  fairly  small  amplitude  that  are  long  in  comparison  with 
the  total  fluid  depth  (e.g.,  see  Benney,  1966).  However,  this  meaning  of 
"long"  is  overly  restrictive  because  it  excludes  internal  waves  whose 
wavelengths  may  be  comparable  or  even  less  than  the  total  fluid  depth,  but 
which  are  much  longer  than  the  thickness  of  an  appropriate  thin  layer 
defined  by  the  background  density  distribution.  For  example,  Osborne  and 
Burch  (1980)  have  observed  internal  waves  in  the  Andaman  Sea  that  seem  to 
behave  like  KdV  solitons,  even  though  their  observed  wavelengths  are  only 
comparable  to  the  total  water  depth. 

As  we  will  discuss  in  Section  2,  the  derivation  of  (2)  permits 
wavelengths  comparable  to  the  total  fluid  depth,  provided  only  that  they  are 
much  longer  than  the  thickness  of  an  appropriate  thin  layer.  Thus,  there  is 
a  sense  in  which  (2)  generalizes  (1),  and  one  might  expect  (2)  to  be  at 
least  as  accurate  as  (1)  in  predicting  experimental  data. 

We  find  the  opposite  to  be  true:  for  our  data  set,  the  predictions  of 
(2)  are  always  less  accurate  than  those  of  (1).  The  finite-depth  equation 
accurately  predicts  the  data  only  Insofar  as  it  agrees  with  the  KdV 
equation.  This  conclusion  Is  based  on  our  limited  set  of  data,  but  Koop  and 
Butler  (1981)  have  reached  substantially  the  same  conclusion  on  the  basis  of 
independent  experiments. 

As  we  will  show,  the  resolution  of  this  paradox  is  as  follows.  Both 

(1)  and  (2)  are  derived  from  Euler's  equations  of  motion  by  comparable 
asymptotic  expansions.  In  each  case,  the  solution  of  the  equation,  (1)  or 

(2) ,  provides  the  dominant  term  in  the  asymptotic  expansion.  However,  In 
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terms  of  allowable  wave  amplitudes,  the  range  of  validity  of  (2)  is  rather 
small;  in  fact.  It  Is  much  smaller  than  that  of  (1).  This  limited  range  is 
found  by  carrying  the  expansion  that  leads  to  (2)  to  the  next  order,  and 
comparing  the  second  order  theory  with  experimental  data  (cf.,  Joseph  and 
Adams,  1981). 

Thus,  under  somewhat  different  conditions,  both  (1)  and  (2)  predict  the 
slow  evolution  of  long  internal  waves  of  small  amplitude  as  they  travel  in 
one  direction.  However,  waves  that  are  long  enough  to  satisfy  the 
requirements  for  (2)  are  not  necessarily  long  enough  for  (1).  On  the  other 
hand,  waves  that  have  amplitudes  small  enough  to  satisfy  the  requirements 
for  (1)  may  lie  outside  the  range  of  validity  of  (2). 


2.  Derivation  of  the  Equations 


We  consider  a  two-fluid  configuration,  in  which  a  layer  of  lighter 
fluid  over  lies  a  layer  of  heavier  fluid,  resting  on  a  horizontal 
Impermeable  bed  in  a  constant  gravitational  field  (see  Figure  2).  This  is 
the  simplest  configuration  that  supports  internal  gravity  waves,  and  it  is 
adequate  to  model  the  waves  observed  in  our  experiments,  in  those  of  Koop 
and  Butler  (1981)  and  in  those  of  Osborne  and  Burch  (1980).  It  excludes 
higher  vertical  modes,  including  those  studied  by  Kubota,  Ko  and  Dobbs 
(1978).  Therefore,  our  results  cannot  be  compared  directly  to  theirs, 
although  the  derivations  themselves  may  be  compared. 

The  problem  of  finding  the  two-dimensional,  infinitesimal,  irrotational 
disturbances  admitted  by  two  stably  stratified  layers  of  incompressible 
fluid  in  a  constant  gravitational  field  was  discussed  by  Lamb  (1932,  Section 
231).  The  velocity  potentials  in  the  upper  and  lower  fluids  may  be  written 


$2  -  (A  sinh  kz  +  B  coshkz)  exp{ik[x  -  c(k)t]}  , 


~  D  sinhk(z+h2)  exp{1k[x  -  c(k)t]}  . 


Here  c(k)  must  satisfy  the  linear  dispersion  relation, 


♦  (1  -  A)  TxT2D 


•  M  ■" 


♦  T2]  ♦  ATjT; 


0  ,  (4a) 
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where 


A  -  ,  (4b) 

P2 

Ti  «  tanh  k  h^  ,  i  ■  1,2.  (4c) 

In  our  experiments  a  «  0.05;  therefore,  we  will  use  the  Boussinesq 
approximation  (a  ♦  0,  but  ga  finite)  to  simplify  results.  Some 

generalizations  to  arbitrary  a(0  <  a  <  1)  are  given  below  and  in  the 

following  paper  (II).  For  small  a,  the  roots  of  (4a)  are 

c$(k)  *  ^  tanh  k(hj  +  h2)  +  0(a)  (surface  waves)  (5a) 

c?(k)  *  -JaTiT2,  ,  [i  +  0(a)]  (internal  waves)  (5b) 

kdi  +  T2) 

In  the  Boussinesq  limit,  these  speeds  always  are  distinct  for  a  given  k. 
Henceforth,  we  will  let  a  ♦  0  (and  g  ♦  •)  and  retain  only  dominant  terms. 

The  major  difference  between  (1)  and  (2)  are  due  to  different 

approximations  of  (5b).  The  KdV  limit  is  obtained  by  letting 

h(hj  +  h£)  ♦  0,  so  that 

Cs(0)  -  g(hi  +  hz)  ,  (6a) 
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< 


< 


9*  hi  h2 
+  h2 


(6b) 


The  dispersive  term  in  (1)  corresponds  to  the  first  correction  to  (6b)  for 
small  k(hj  ♦  h2).  On  the  other  hand,  the  finite-depth  limit  amounts  to 
h2/hi  +  0.  khi  finite,  so  that 

c?(0)  *  9Ah2  £  c 1  (7) 

The  dispersive  term  in  (2)  corresponds  to  the  first  correction  to  (7)  for 
small  kh2*  Note  that  the  speeds  of  the  long  infinitesimal  waves  in  (6b;  and 
(7)  differ,  unless  we  also  require  f^/hj  ♦  0  in  (6b). 


2.1.  The  KdV  Equation 

Aspects  of  KdV  theory  for  long  internal  waves  have  been  discussed  by 
Keulegan  (1953),  Long  (1956),  Peters  and  Stoker  (1960),  Benjamin  (1966), 
Benney  (1966)  and  others  (cf..  Miles,  1979,  and  the  references  cited 
therein).  We  now  briefly  outline  the  derivation  of  (1)  for  the  two  layer 
configuration  shown  in  Figure  2.  The  basic  assumptions  are: 

1)  The  waves  are  long  In  comparison  with  the  total  fluid  depth: 

k2  (hA  +  h2)2  «  1  , 
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where  k-1  represents  a  characteristic  horizontal  wavelength. 

11)  The  waves  are  small,  so  that  If  *n  denotes  a  characteristic  wave 
amplitude,  then 


n/(hx  +  h2)  «  1  . 


ill)  The  two  effects  are  In  approximate  balance: 

e  *  nAhj  +  h2)  «  0Ck2(h2  +  h2)2]  «  1  . 


iv)  Viscous  effects  are  weaker  than  either  of  these. 

(In  addition,  we  assume  throughout  this  series  that  the  motion  is 
two-dimensional,  and  that  the  fluid  is  incompressible.) 

It  is  consistent  with  these  assumptions  to  define  dimensionless  (*) 
variables  as  follows 


z* 


z 

h  j  +  h2 


x* 


C1'2X 

hx  +  h2  * 


We  also  Introduce  a  slow  time  variable, 


t*  ■  c  t*  . 

The  Interface  Is  defined  by 


(9) 
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n  «  t(h^  +  h2)  n*(x*.  t*.  t*;  e)  ; 


if  the  upper  surface  is  free,  it  is  defined  by 


(10a) 


C  -  hx  +  c(hx  +  h2)  c*(x*.  t*.  t*;  c) 


(10b) 


There  is  a  velocity  potential  in  each  layer.  In  the  lower  fluid,  because  of 
Laplace's  equation  and  the  boundary  condition  at  z  *  -hg,  the  potential  has 
the  formal  expansion. 


*2  •  *o(x*.  **»  e)  -  Y  \z*  4  h 


h  2  V  32»q 
i  ♦  b2j  a(x*)2 


+  0(c4)  . 


There  is  a  corresponding  expansion  in  the  upper  fluid.  At  the  interface, 
the  normal  velocity  and  the  presssure  must  be  continuous,  while  the  pressure 
must  vanish  at  the  free  surface.  In  addition,  both  the  free  surface  and  the 
interface  satisfy  kinematic  conditions  (Dc/Dt  ■  Dn/Dt  •  0).  Finally,  all 
motion  should  vanish  as  |x|  ♦  •.  In  order  to  satisfy  these  conditions 
order- by-order  in  c,  it  is  necessary  to  expand  n*,  £*,  and  the  velocity 
potentials.  Thus,  for  example 


r»*  *  *  cn2* 
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At  the  leading  order  (infinitely  long  waves  of  infinitesimal 
amplitude),  the  equations  are  hyperbolic  and  linear.  An  initial  disturbance 


Is  decomposed  into  four  wave  modes: 


where  [cj]  are  the  four  roots  in  (6).  Because  these  speeds  are  so 
different,  localized  surface  and  internal  waves  quickly  separate  in  space. 

Weak  nonlinear  interactions  and  weak  dispersive  effects  appear  on  the 
next  scale  [t*  *  0(1)],  when  the  expansion  is  carried  to  the  next  order  and 
secular  terms  are  eliminated.  There  are  no  nonlinear  interactions  between 
modes,  because  they  interact  with  each  other  for  too  short  a  time.  However, 
each  mode  interacts  with  itself  for  a  long  time,  and  each  of  the  four  waves 
in  (11)  evolves  according  to  its  own  KdV  equation.  Omitting  details  of  the 
analysis,  the  dimensional  equation  for  the  internal  wave  mode  is 


c  i  in  +  in  +  1  A  _i\  in 
0  at  ax  2  ^h2  hlJ  ax 


+ 


l 

6 


0  , 


(12) 


where  co  is  given  by  (6b).  [Recall  that  (12)  is  valid  for  a  «  1.  The 
generalization  to  arbitrary  a  is  given  in  the  following  paper  (II)  if  the 
upper  surface  is  free,  and  by  Djordjevic  and  Redekopp  (1978)  if  the  upper 
surface  is  rigid.]  To  reduce  (12)  to  (1),  let 


< 


X 


X  -  c0t 

(PPP7-  T 


Then  f(x.t)  satisfies 


f,  ♦  6ffx  ♦  fm  -  0  . 


We  now  state  some  of  the  consequences  of  KdV  theory  that  may  be  tested 
experimentally.  More  details  may  be  found  in  Segur  (1973),  Hammack  and 
Segur  (1974,  1978),  and  elsewhere. 

i)  The  soliton  solution  for  KdV  is 


f(x.T)  -  2k2  sech2  (<c(x  -  4ie2  t  -  Xq)}  (14a) 

where  <c,xo  are  arbitrary  constants.  In  dimensional  terms, 

n  ■  n  sech2  (p(x  -  pvt  -  x0)}  (14b) 

where 


3  (h1  -  h2) ~ 

4  hfh2 


1  hi  -  h 

2  hj  h 


(14c) 
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These  results  were  given  first  by  Keulegan  (1953),  except  for  a 
misprint.  Because  f  >  0,  it  follows  that  an  internal  soli ton  always 
thickens  the  thin  layer.  Thus,  the  soliton  raises  the  interface  if 
hj  >  h2,  and  lowers  it  if  hi  <  h2* 

ii)  In  our  experiments,  n  >  0  initially.  It  follows  that  the  internal 
wave  evolves  into  soli  tons  when  hi  >  h2,  but  not  when  hi  <  h2* 
ill)  The  nonlinear  term  in  (12)  vanishes  if  hi  »  h2»  We  will  examine  this 
special  case  in  detail  in  part  III  of  this  series.  Here  we  assume 

hi  i  h2* 

iv)  Arbitrary  initial  data  that  are  smooth  and  localized  will  evolve  into 
N  solitons,  ordered  by  amplitude,  followed  by  a  dispersion  oscillatory 
tail.  The  number  of  solitons  that  emerge  from  f(x.O)  is  the  mxnber  of 
zeros  of 


2 

+  T(x»0)  v*  ji  ♦  1  is  X  ♦  **  (1®) 

dx 


2.2  The  Finite  Depth  Equation 

Kubota,  Ko  and  Dobbs  (1978)  derived  (2)  for  cases  in  which  the  fluid  is 
confined  between  two  rigid  walls  and  the  background  density  distribution  is 
continuous.  The  derivation  we  present  next  differs  slightly  from  theirs 
because  our  background  density  distribution  Is  discontinuous.  A  more 
Important  difference,  however.  Is  that  our  choice  of  small  parameters 
(c  «  1)  differs  from  theirs.  Both  derivations  lead  to  (2),  but  the  scaling 
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of  the  physical  variables  is  somewhat  different.  For  simplicity,  we  will 
use  the  Boussinesq  approximation,  and  also  replace  the  free  surface  with  a 
rigid  lid. 

The  assumptions  underlying  (2)  are  that: 

i)  there  is  a  thin  (lower)  layer, 

e  »  h2/hi  «  1; 

ii)  the  characteristic  horizontal  wavelength  is  comparable  to  the  depth  of 
the  thick  layer. 


khl  «  0(1); 


iii)  wave  amplitudes  are  small, 


n/h2  «  1; 


iv)  these  two  effects  balance. 


n/h2  »  0(c); 


v)  viscous  effects  may  be  neglected. 

Note  that  the  assumption  of  long  waves  (kh2  «  1)  is  implied  by  (i)  and 
(ii).  Note  further  that  assumption  (iv)  is  not  consistent  with 
corresponding  assumption  for  KdV.  Consequently,  we  may  expect  the  KdV  limit 
of  (2)  to  be  somewhat  singular. 

The  following  scaling  is  consistent  with  the  assumptions  listed  above. 
In  the  lower  layer,  dimensionless  (*}  variables  are: 
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x*  ■  x/hj  ,  z*  »  z/h2  ■  z/ehx  ,  t*  ■  ct/hx  , 


et*  .  (16a) 


the  wave  speed,  c.  Is  to  be  determined.  The  interface  is  defined  by 


n  ■  eh2Cni  +  en-*  +  e2^]  ♦  0(e4)  . 


(16b) 


After  satisfying  Laplace's  equation  in  the  lower  fluid  and  the  boundary 
condition  at  z  «  -h2,  one  finds  that  the  velocity  in  the  lower  field  at  the 
interface  may  be  represented  by 

u  ■  eTCu^x^.t*;  t,c)  +  eu2]  +  0(e3)  , 


-  e2c 


[3U;  au 

ax*  +  c  ax 


i*  ♦  0(c4)  . 


Because  no  further  terms  will  be  needed  in  this  expansion  in  order  to  derive 
(2),  it  will  follow  that  the  entire  dispersive  effect  in  (2)  is  due  to  the 
upper  fluid.  It  will  also  turn  out  that  the  entire  nonlinear  effect  is  due 
to  the  lower  fluid. 

The  velocity  potential  in  the  upper  layer  may  be  written  as 


1  7  cosh  k(h,  -  z)  ,  #J1  .  .. 

♦>  ■  T,.{  — Slnhk'hr  «*Pdkx)<lk  . 
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It  is  consistent  with  assumption  (ii)  to  define 


m  »  khx  ,  x*  ■  x/hj  , 


to  rewrite  4m  in  terms  of  these  variables.  Then  if  we  expand  if  as 


¥  •  h|  c‘{A1  +  cA2  +  •••)  * 


the  velocites  in  the  upper  fluid  at  the  interface  take  the  form 


U  „  f  (m  cothmJCAi  +  cA2  +  ...)  exp(imx*)dm 


w  ..  -  f  m(AA  +  cA2  +  ...)  exp(imx*) dm  , 


There  are  two  kinematic  conditions  at  the  interface: 


in  +  u  *  w  , 
at  ax 


(2la) 


(U-u)  •  W  -  w  • 


(21b) 
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Continuity  of  pressure  at  the  Interface  may  be  written  in  terms  of 
Bernoulli's  law, 

ill  *  M  *  *‘l  -  ito  -  {-u2  -  so  •  constant  ,  (21c) 

at  2  at  2 

once  we  ha/«  made  the  Boussinesq  approximation.  Because  (21c)  is  valid 
along  the  interface,  its  tangential  derivative  vanishes  there.  Substituting 
(17)  and  (20)  into  these  conditions,  one  finds  that  non-trivial  solutions 
exist  at  leading  order  if 


c2  «  gah2  ,  (22) 

as  anticipated  by  (7).  [A  major  difference  between  this  derivation  and  one 
following  Kubota,  Ko  and  Oobbs  (1978)  is  that  their  leading  order  wave  speed 
satisfies 


■£■2  ■  ga  h2(l  -  e)  . 


Thus,  they  retain  a  term  at  leading  order  that  we  regard  as  a  higher  order 
effect.  This  difference  in  ordering  persists  at  higher  orders  in  the 
expansion.]  In  the  present  derivation,  the  solution  of  the  leading  order 
equations  are 

ni(x*.t*,x)  »  f(r,t)  ♦  g(i.t) 
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Ui(x*,t*,i) 


f(r,t)  -  g(i,T) 


(23) 


Ai(x*,t*,T)  ■  if(m,-r)  exp(-imt*)  -  ig(m,T)  exp(imt*) 


where 


r  ■  x*  •  t*  ,  t  ■  x*  ♦  t*  » 


(*)  denotes  the  Fourier  transform  (assumed  to  exist)  and  we  have  used  the 
boundary  condition  that  all  motion  ceases  as  |x|  *  •. 

Secular  terms  arise  at  the  next  order  unless  the  right-going  waves 
satisfy 


2~  +  3f-|"-'n:T~  cothm  f(m,i)  exp(imr)dm  «  0  ,  (24a) 

OT  O'  •  0i  J 


and  the  left-going  waves  satisfy  a  similar  equation.  This  nonlinear 
evolution  equation  has  several  representations.  Another  is 


*  3f£  -  1$  /'Othf  (r-y)f(y.T)dy  .  0.  (24b) 


which  may  be  scaled  to  (2).  In  dimensional  variables, 
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/rr\_  I  in  .  in  3  n  3n  Jh  82  f  *  x  -  z  Jz\  .  n 

»  +  w  *  I  hi  »  •  4hl  I?  /  coth2  ~  ihlj  dz  °  • 


(25) 


where  o’  is  defined  by  (22). 

In  this  derivation  we  have  assumed  that  all  functions  vanish  as 

|x|  ♦  •.  Periodic  boundary  conditions  are  more  natural  in  many  problems, 
such  as  when  the  waves  are  generated  by  the  periodic  motion  of  the  tides. 
In  this  case  it  is  necessary  only  to  replace  the  Fourier  integrals  with 
Fourier  sums.  The  analogue  of  (24)  on  (-*,»)  with  periodic  boundary 

conditions  is 

2  ■—  +  3f~-  -  n  coth  n  fn(t)  exp(inr)  dn  «  0  ,  (26a) 

ar  ar 

where  [fn]  are  the  coefficients  of  the  Fourier  series  representation  of  f. 
The  sum  in  (26a)  can  be  written  as  a  convolution  integral,  with  a  kernel 
involving  elliptic  functions.  The  derivation  of  (26a)  also  requires 

« 

f°{  t)  ■  J  f(r,x)  dr  ■  0  ,  (26b) 


which  amounts  to  a  normalization.  The  analysis  of  (26)  Is  quite  similar  to 
that  of  (24),  as  discussed  by  Ablowitz,  Satsuma,  Fokas,  and  Segur  (1981). 
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Based  on  the  work  of  Joseph  and  Egrl  (1978),  Chen  and  Lee  (1979), 
Satsuma,  Ablowitz  and  Kodama  (1979),  and  Kodama,  Satsuma  and  Ablowitz 
(1981),  we  may  assert  that  (24)  is  completely  integrate,  and  that  Its 
solutions  are  qualitatively  similar  to  those  of  the  KdV  equation.  In 
particular,  a  soliton  solution  of  (24)  is 

3  .  _  x  sin  X _ _  (27a) 

2  ‘  *  cos  x  +  coshx[r  ♦  (t/2)  x  cot  x  +  r0] 

where  (x,ro)  are  arbitrary  parameters  except  that  0  <  x  <  *•  In  dimensional 
terms. 


JL  ,  (2/3)(h2/h1)  x  sin  x 

h2  *  cos  x  *  cosh  e 


(27b) 


,  .  x  («  -  x°>.  .  [1  .  i  x  cot  0  (27c) 

hj  hj  Z  hj 

This  reduces  to  a  KdV  soliton  in  the  limit  X  ♦  0,  h2/hl  fixed. 

In  order  to  carry  this  expansion  to  higher  order,  it  is  convenient  to 
consider  only  right-going  waves,  so  that  g(4,t)  =  0  in  (23).  Then  at  second 


order. 


U  ■  7  <2  *  I^TCf]  . 


2^  f  A2  exp(1mx*)  dm  *  -f2  -  |f2  -  . 


At  second  order  f2(r,-r)  is  free,  but  secular  terms  arise  at  third  order 
unless 


««  •  *  H(f3> 


l  iff 
3  ar*  ’ 


where  f  is  a  solution  of  (24),  and 


L[v]  -  2  -fi  +  3  (fv)  ♦  T[v]  . 


Because  L(v)  is  the  linearization  of  (24),  one  finds  easily  that 


LE**)  ■  (.[“]  ■  o  . 

3r  3t 


(30a) 


(30b) 


Ufl  -  |£(f*>  • 


(30c) 
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We  now  seek  specifically  the  second-order  correction  to  a  soliton  (27) 
so  that 


if 

3t 


-  (Cl  +  ec2)  , 


i£& 

3t 


and  c2  is  to  be  determined.  Moreover, 


£t[f]  ♦!  fl  -  2c  i  f  -  0 


and 


«  X2f  +  9Clf2  .  If3 
In  this  case,  (29)  reduces  to 

l[f2]  »  jp  {(2c2  -  3c2  -  ~)  f  -  3c j  f2  ♦  \  f3}  .  (31) 

From  (30a),  this  equation  has  no  unique  solution.  From  (30b),  we  must 
choose 


2c2  •  3cf  ♦  \2/3 


(32) 


in  order  to  avoid  unbounded  growth  in  f2.  Thus  the  dimensional  speed  of  a 


soli ton  to  this  order  Is 


C  ~  c  [1  -  |  ex  cot  X  +  (3/4  cot2  X  +  1/3)]  ,  (33) 

where  c  «  h2/hl,  and  X  Is  the  parameter  for  the  soliton. 

In  the  small  amplitude  limit,  X  ♦  0  and  (33)  reduces  to  the  expansion 
of  (6b)  for  small  e.  At  the  other  extreme,  if  x  it  then  to  this  order  of 
approximation. 


C 


T  -  X/ 


-1/2 


Clearly 


X  <  It(l-e)  (34) 

is  necessary  for  the  validity  of  (33)  as  an  asymptotic  expansion.  It 
probably  is  not  sufficient.  Joseph  end  Adams  (1981)  suggest  that  their 
expansion  is  valid  only  if  C  <  1.4c*,  where  c*  corresponds  to  Co.  defined  by 
(6b).  In  this  problem,  their  cutoff  corresponds  to  about 

X  <  u(l  -  2c)  . 

However,  they  give  no  clear  criterion  for  this  choice. 


We  now  return  to  (31).  Let  v(r)  denote  a  particular  solution  of 


L[v]  -  (f3)  . 

or 


(35) 


We  have  found  no  analytic  expression  for  v(r),  but  a  numerical  procedure  to 
solve  (35)  approximately  is  outlined  in  an  appendix.  A  solution  of  (31), 
corresponding  to  f[r  -  (c2  +  c  c2)t;a]  is 


*2 


-  2c,  f  *  \  ' 


(36) 


or 


n  ~  ch2  [(1  -  2e  Cj)f  +  ^  ev]  . 

This  is  the  expression  we  will  compare  with  our  data  in  Section  3. 

In  the  experiments  of  Koop  and  Butler  (1981),  the  dimensionless  density 
difference  a  «  0.367,  and  the  Boussinesq  approximation  is  invalid.  However, 
the  analysis  for  arbitrary  a  follows  nearly  identical  lines,  and  we  may 
simply  state  the  main  results.  With  the  thin  layer  on  the  bottom  and  a  free 
surface  on  top,  the  leading  order  wave  speed  is  still  given  by  (22).  The 
evolution  equation  becomes 


2£  *  «£  ♦  ■  ° 


instead  of  (24).  For  fixed  X,  the  relation  between  a  soliton  with  &  4  0(f, 
cx)  and  the  earlier  results  (f,  cj  is 

CiU;  L)  ■  (1  -  a)  c i ( X )  ,  ?  ■  (1  -  a)  f(r  -  Cit;  X)  .  (38) 

The  evolution  equation  for  f2  is  somewhat  complicated,  but  if  f  is  a 
soliton  and  ?2  is  a  permanent  wave  traveling  with  the  same  speed,  then 


CCf2]  5  -£l-2c,  f2*  3ffj  ♦  (l-4)-^TC?2]) 


3c  22  *  7  ~  1 7 a  +  13  a2  3a  3  Tr*i  \ 

T  2(1  -  A)2  ■  2  3r  L  J  » 


(1  -A); 


The  result  is 
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h 


•  -  2c ! 


(1  *  A  -  A2) 

(l-*)2 


? 


3  A 

2(1  -  A) 


?2  ♦ 


[2(1  -  A)*  +  |]  v  , 


where  (f,  cx)  satisfy  (32),  and  v  satisfies  (31).  Let  (n)  denote  the 
maximum  displacement  of  the  interface.  Then  it  follows  from  (40)  that 

hjTi’  (2/3)x  sinx  f  1+A-A2  eAXSinxl 

Tl  “  A ) h |  "  1  +  COS  X  |1+eXCOtX  (1-A)2  ”  1  +  COS  X  J 


+ 


e  2(1-  A)  + 


2(1  -  A) J  max 


0(e2) 


(41) 


Koop  and  Butler  (1981)  define  an  integral  length  scale  by 


(42) 


It  follows  from  the  results  above  that 


j|  -  ^  (1  -  A)  +|  eX2  COt  X  +  c  (2(1  -  A)2  +  |  j  I  +  0(e2)  ,  (43) 


where 


r 

1 


V(r)  dr 


i 
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and  must  be  found  nunerlcally. 
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3.0  Comparison  with  Experiments 


A  series  of  experiments  was  performed  in  which  long  internal  waves  were 
generated  by  the  vertical  uplift  of  a  rectangular  piston  at  one  end  of  a 
wave  tank,  as  shown  in  Figure  1.  The  tank  was  30  m  long,  60  cm  deep,  and 
39.4  cm  wide.  The  piston  was  61  cm  long  and  spanned  the  tank  width  (see 
Hammack,  1972,  for  details).  The  end  wall  of  the  tank  near  the  piston  acted 
as  a  plane  of  symmetry  (in  the  linear,  inviscid  theory)  so  that  the 
effective  elevation  of  the  interface  was  twice  as  long  (122  cm)  as  the 
piston,  with  a  correspondingly  decreased  amplitude.  In  these  experiments, 
hi  «  45  cm,  h2  *  5  cm,  and  A  ■  0.048.  The  thickness  of  the  pycnocline 
gradually  increased  from  about  1  cm  to  2  cm  during  the  experiments.  The 
total  displacement  of  the  piston  was  2  cm. 

Impulsive  upward  motion  of  the  piston  from  the  bed  of  the  tank 
generated  both  suface  and  internal  waves.  [See  Hammack,  1981,  for  an 
analysis  of  the  generation  process.]  The  surface  wave  moved  so  much  faster 
than  the  internal  wave  that  there  was  effectively  no  interaction.  A  large 
vertical  plate  that  spanned  the  tank  width  was  carefully  lowered  into  the 
water  at  a  point  18.8  meters  downstream  of  the  piston  after  the  surface  wave 
had  passed  that  point.  This  plate  prevented  the  surface  wave  from 
re-entering  the  test  section,  where  the  much  slower  Internal  wave  was  still 
propagating.  Moreover,  It  effectively  lengthened  the  test  section,  since 
the  internal  wave  eventually  reflected  off  the  plate  and  propagated  back 
through  the  test  section.  * 
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Internal  wave  measurements  were  measured  using  a  laser-optics-detector 
system  described  by  Hammack  (1981).  Basically  the  light  beam  of  a  laser  was 
converted  to  a  uniform-width  light  sheet  that  was  directly  across  the 
glass-walled  tank.  With  the  quiescent  interface  intercepting  the  light 
sheet  and  the  brine  dyed  blue,  subsequent  motions  of  the  interface  varied 
the  amount  of  light  traversing  the  tank  cross-section.  The  transmitted 
light  was  focused  onto  a  photo-detector  which  produced  an  output  voltage 
proportional  to  the  amount  of  transmitted  light.  The  system  was  calibrated 
before  each  experiment  to  determine  the  (nonlinear)  correlation  between 
output  voltage  and  vertical  displacement  of  the  interface.  Only  one  gage 
was  available;  hence,  the  same  experiment  was  repeated  in  order  to  measure 
time  histories  of  the  interfacial  displacement  at  various  positions  along 
the  tank.  Both  the  repeatability  of  the  piston  motion  and  the  quiescent 
conditions  prior  to  an  experiment  were  carefully  monitored  and  assured  for 
all  reported  data. 

Figure  2  shows  the  evolution  of  the  internal  wave,  as  measurued  at 
seven  downstream  locations  in  the  tank.  The  coordinate  system  in  Figure  2 
moves  with  the  linearized  wave  speed,  co,  and  reverses  the  sense  of 
direction  of  the  wave.  Thus,  the  front  of  these  waves  is  to  the  left  in 
each  record.  Moreover,  we  have  omitted  the  weak  train  of  oscillatory  waves 
that  follows  the  dominant  waves  that  are  shown.  The  first  five  measurements 
were  recorded  ahead  of  the  plate  at  x/(h1h2)i'2  ■  125. 

For  this  experiment  we  may  take  k"1  ■  122  cm,  and  "n  ■  1  cm, 
corresponding  to  the  dominant  length  scale  and  maximum  amplitude  of  the 
Initial  wave.  Hence,  k2(h;  +  h2)2  ~  0.16,  n/(hj  ♦  h2)  -  0.02,  and  we  may 
apply  KdV  theory  at  least  provisionally.  Using  either  the  first  or  second 
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wave  record  as  the  potential  in  (15),  one  finds  that  according  to  KdV 
theory,  this  wave  should  evolve  into  two  solitons,  followed  by  a  weak 
oscillatory  wave  train.  Certainly  this  prediction  is  in  qualitative 
agreement  with  the  wave  records  shown  in  Figure  2. 

A  more  rigorous  test  of  KdV  theory  is  shown  in  Figure  3,  where  we  have 
plotted  the  lead  wave  from  each  of  the  last  five  records  on  a  single  graph. 
According  to  (14),  the  data  from  all  of  these  wave  records  should  fall  on  a 
single  curve  once  the  peak  amplitude  at  each  station  is  known.  The 

agreement  in  Figure  3  between  the  predicted  wave  shape  and  all  of  these  data 
is  so  good  that  we  conclude  that  these  are  (locally)  KdV  solitons.  The 
major  discrepancies  between  the  predicted  and  observed  wave  shapes  occur: 
(i)  at  the  rear  (i.e.,  right)  of  each  wave,  where  the  influence  of  the 

trailing  soliton  becomes  important;  and  (ii)  in  the  data  at 

x/(hlh2)l/2  •  151,  just  downstream  of  the  reflecting  plate,  where  the 

incident  and  reflected  waves  from  this  plate  apparently  are  still 

superposed. 

The  amplitude  of  the  soliton  is  slowly  decaying  as  it  propagates  down 
the  tank,  presumably  due  to  viscosity,  but  this  decay  is  sufficiently  slow 
that  the  wave  continually  readjusts  its  shape  as  it  decays  so  that  it  is 
locally  a  KdV  soliton.  On  the  basis  of  a  much  larger  set  of  experiments, 
Hammack  and  Segur  (1974)  found  that  this  same  description  (locally  KdV 
solitons,  slowly  attenuated  by  viscosity)  also  applies  to  the  corresponding 
long  surface  waves.  See  the  following  paper  (II)  for  a  discussion  of  the 
viscous  decay  of  these  waves. 


31 


A 


Based  on  the  remarkably  good  agreement  in  Figure  3,  and  fact  that 
k2(h2  +  h2)2  ~  0.16  initially  for  these  data,  it  is  tempting  to  conclude 
that  the  KdV  equation  has  a  fairly  large  range  of  validity.  This  conclusion 
may  well  be  correct,  but  it  is  not  necessarily  implied  by  Figure  3.  While 
k2(hi  +  h2)2  ~  0.16  may  describe  the  initial  data,  it  does  not  represent  the 
evolving  solitons,  whose  characteristic  wavelengths  are  much  longer.  The 
stronger  conclusion  applies  only  if  the  equation  correctly  predicts  from  the 
initial  data  the  wavelengths  of  the  solitons  that  emerge. 

KdV  theory  predicts  that  all  solitons  should  move  somewhat  faster  than 
Co,  the  linear  long  wave  speed,  and  therefore  should  move  slowly  to  the  left 
in  Figure  2.  The  larger,  faster  soliton  does  move  somewhat  to  the  left 
until  the  last  station,  but  the  smaller  soliton  clearly  is  moving  slowly  to 
the  right.  This  behavior  also  was  observed  for  surface  solitons  (Hammack 
and  Segur,  1974).  It  seems  to  be  a  viscous  effect  that  all  solitons  travel 
slower  than  the  prediction  of  KdV  theory,  and  often  slower  than  Co. 

Next,  we  consider  the  finite-depth  equation,  (24).  There  seems  to  be 
no  simple  way  to  compare  all  of  the  data  on  one  graph  for  (24),  but  in 
Figure  4  we  have  replotted  the  lead  waves  measured  at  two  representative 
points,  x/(hih2)l/2  ■  60  and  191.  This  figure  also  shows  the  wave  shapes 
predicted  by  KdV(14),  by  the  first-order  finite-depth  theory  (27),  and  by 
the  second-order  finite-depth  theory  (36).  In  all  cases,  the  free  parameter 
available  in  the  theory  [e.g.,  x  in  (27)]  was  chosen  to  match  the  peak 
amplitude  of  the  measured  wave.  It  is  apparent  that  while  all  of  these 
theoretical  predictions  are  reasonably  accurate,  first-order  finite-depth 
theory  Is  noticeably  less  accurate  than  KdV  theory.  The  discrepancy  Is 
apparent  even  at  x/(hih2)1/2  •  191,  where  h2/hl  -  1/9,  lj/h2  -  0.094  and  one 
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might  expect  (27)  to  be  quite  accurate. 

Figure  4  supports  the  claim  that  the  range  of  validity  of  the 
asymptotic  expansion  that  leads  to  the  finite-depth  equation  is  small.  In 
Figure  4a,  where  the  peak  wave  amplitude  is  rather  small,  the  second-order 
finite-depth  theory  predicts  the  measured  data  about  as  well  as  does  KdV 
theory.  The  peak  amplitude  in  Figure  4b  is  about  twice  that  in  Figure  4a, 
and  here  even  the  second-order  theory  is  beginnning  to  fail.  Presumably, 
third-order  corrections  now  have  become  important. 

As  an  independent  test  of  the  hypothesis  that  the  range  of  validity  of 
(24)  is  rather  small  in  a  practical  sense,  we  analyze  next  the  experimental 
results  of  Koop  and  Butler  (1981).  Figure  5  shows  the  relation  between 
soliton  amplitude  (41)  and  integral  length  scale  (43),  corresponding  to 
Figures  10  and  12  of  Koop  and  Butler.  The  first-order  curve  is  obtained  by 
letting  e  ♦  0  in  (41)  and  (43).  Higher-order  corrections  depend  on 
<(■  h2/hl)  and  on  A  ■  (1  -  pi/p2),  and  we  have  plotted  in  Figures  5a,  5b  two 
second-order  curves,  corresponding  to  the  two  configurations  of  Koop  and 
Butler.  Note  that  both  second-order  curves  terminate  at  finite  values  of 
Tihi/(1  -  A)h£.  This  termination  occurs  because  hhl/(l  -  A)h£  has  a  maximum 
in  (41)  as  a  function  of  A,  signalling  a  breakdown  of  the  asymptotic  series. 
This  breakdown  occurs  earlier  than  that  in  (34). 

A  comparison  of  the  first-order  theory,  second-order  theories,  and  the 
data  in  Figure  5  reveals  the  following  facts. 

1)  In  Figure  5a  (e  ■  0.197,  corresponding  to  Figure  10  of  Koop  and 
Butler)  there  are  no  data  with  amplitudes  small  enough  that  the  first  and 
second-order  theories  coincide.  In  this  sense,  the  first-order  theory  by 
Itself  must  be  considered  Inadequate  to  represent  this  set  of  data. 
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2)  For  e  ■  0.197,  the  second-order  theory  predicts  the  data  quite  well 
within  Its  range  of  validity.  However,  a  significant  portion  of  the  data 
lies  outside  the  range  of  validity  of  the  theory. 

3)  In  Figure  5b  (e  *  0.029)  the  flrst-and  second-order  theories  are  in 
close  agreement  over  the  entire  range  of  validity  of  the  theory.  Both 
predict  wave  amplitudes  somewhat  larger  than  those  observed,  especially  for 
the  longer  (and  therefore,  smaller)  waves.  This  effect  is  less  pronounced 
but  also  evident  In  Figure  5a.  Koop  and  Butler  noted  that  even  the  KdV 
equation  has  this  problem  for  very  long  waves;  they  argued  that  It  is  a 
viscous  effect. 

4)  For  c  *  0.029,  most  of  the  data  lies  beyond  the  range  of  validity 
of  the  theory.  There  is  no  obvious  physical  mechanism,  such  as  wave 
breaking,  that  precipitates  the  breakdown  of  the  theory. 

Finally,  let  us  summarize  our  major  conclusions. 

1)  Both  (1)  and  (2)  are  valid  asymptotic  equations  that  govern  the 
slow  evolution  of  long  internal  waves  of  small  amplitude  that  propagate  in 
one  direction  in  an  inviscid  fluid.  However,  the  meaning  of  "long"  is 
different  in  the  two  theories,  which  describe  perturbations  to  two  different 
wave  equations,  with  two  different  speeds. 

2)  The  KdV  equation  (1)  seems  to  have  a  relatively  large  range  of 
validity,  and  has  practical  predictive  value  even  when  its  requirements  are 
satisfied  only  marginally. 

3)  The  finite  depth  equation  (2)  by  itself  has  such  a  small  range  of 
validity  that  it  has  been  difficult  even  to  find  it  experimentally. 
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4)  The  asymptotic  expansion  that  generates  (2)  has  a  larger  range  of 
validity,  but  it  is  also  quite  limited.  In  particular,  there  is  a  range  of 
wave  amplitudes  for  which  the  long  waves  in  question  are  too  small  to  break, 
but  too  large  to  be  predicted  by  the  finite-depth  theory,  of  any  order. 
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Figure  Captions 


1.  Piston  moves  up  Into  a  two-layer  system,  generating  waves  both  at  the 
free  surface  (;)  and  at  the  Interface  (n)» 

2.  Internal  waves,  measured  In  time  at  seven  successful  locations, 
hi  «  45  cm,  h2  •  5  cm,  A  ■  0.048.  A  vertical  plate  inserted  during  the 
experiment  at  x/(hih2)l/2  ■  125  reflected  the  wave  train  back  Into  the 
test  section  for  the  last  two  measurements. 

3.  Comparison  of  the  lead  waves  in  the  last  five  measurements  in  Figure  2 

with  the  shape  of  a  KdV  soliton,  according  to  (14). 

C2  «  [3ga(hi  -  h2)/4(hlh2)3]  *  8063. 

4.  Comparison  of  two  of  the  measured  lead  waves  in  Figure  2  with  three 
theoretical  shapes,  each  having  the  same  peak  amplitude  as  that 

measured:  — — ,  KdV  (14);  _ _ _ _  1st  order  finite-depth  (26); 

_ _  2nd  order  finite-depth,  (36);  measured. 

5.  Relation  between  integral  wavelength  (43)  and  peak  amplitude  (41)  for 

internal  solitary  waves  in  a  fluid  of  finite  depth. 

(a)  e  »  h2/hj  »  0.197,  a  ■  1  -  px/p2  *  0.367;  _ _  1st  order  finite 

depth  theory,  with  e  ■  0  in  (41)  and  (43);  — — ,  2nd  order  finite  depth 
theory;  0,  incident  waves  measured  by  Koop  and  Butler  (1981);  t, 

measured  reflected  waves,  (b)  c  .»  0.029,  a  ■  0.367;  - -  1st  order 

finite-depth  theory;  _ - _ _  2nd  order  finite-depth  theory;  a, 

measured  incident  waves,  a,  measured  reflected  waves. 
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Appendix 


For  permanent  localized  waves,.  (35.)  is  equivalent  to 


L  [v]  s  (3f  -  2c)  v  +  -2-  T[v]  -  f3  ,  (A.l) 

l  l  ar 

where  f(r;A)  denotes  a  soli ton  with  speed  Cj(x).  A  Galerkin  procedure  may 
be  used  to  solve  (A.l)  approximately.  Here  (A.l)  is  replaced  by  a  finite 
set  of  equations  of  the  form 

(LjCvD.  ♦„)  •  (f3.  ♦„)  ,  n  •  0.....N  .  (A. 2) 

where  (a,b)  ■  f  ab  dr,  and  Un(r)3  denotes  a  set  of  basis  functions  of  a 
finite-dimensional  space  of  functions  on  -•  <  r  <  ».  We  have  chosen  to  use 

♦n(r)  «  Hen(r)  »exp(-r2/4)// nT  (A.3) 

where  Hen(r)  are  Hermite  polynomials  (notation  as  In  Abramowitz  and  Stegun, 
1964).  These  basis  functions  could  be  generalized  by  Including  a  scaling 
factor  (r  ♦  or),  but  (A.3)  was  adequate  for  our  purposes.  These  functions 
are  orthogonal,  and  their  Gaussian  decay  mimics  the  localized  nature  of  the 
sollton.  Therefore,  If  the  sollton  width  roughly  matches  the  width  of  the 
Gaussian  filter,  then  very  few  basis  functions  are  required  to  represent  a 
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soliton  to  a  high  degree  of  accuracy. 

A  soliton  is  an  even  function  in  Its  argument.  Both  LlCO  and  the 
right-hand  side  of  (A.l)  also  are  even,  and  (A.l)  has  a  solution  that  is 
even.  Therefore,  we  use  only  the  even  basis  functions,  4>2n(r),  in  (A. 2). 
This  choice  excludes  the  solution  of  the  homogeneous  problem  in  (30a),  and 
seems  to  make  (A. 2)  well-posed.  Certainly  it  is  computationally  stable  if 
only  [$2n]  are  used. 

After  expanding 

f(r)  ■  I  f2„  Wr)  ’  ,(r>  ’  |  52n  Wr>  •  <» 

(A. 2)  reduces  to  a  set  of  linear  algebraic  equations  for  [$2n],  the 
coefficients  of  which  equations  being  obtained  by  evaluating  certain 
Integrals.  The  only  delicate  question  is  to  decide  what  accuracy  is 
achieved  with  a  certain  N  in  (A.4).  We  used  three  main  tests  of  nimerical 
accuracy. 

1)  For  a  soliton,  f(r;x),  we  have  that 


Using  (A.4)  In  (A. 5)  yields  a  sequence,  2cx(x»  N),  which  may  be  compared  to 
the  exact  result,  2cx(x)  ■  -xcotx.  With  N  ■  7,  the  error  in  (2ci)  remained 
below  10-2  for  «/10  <  x<  */2,  Increased  to  about  0.2  for  x  •  2«/3,  and 
continued  to  Increase  for  larger  values  of  x.  Regardless  of  how  well  it 
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approximates  (-xcotx),  the  solution  of  (A. 5)  is  the  appropriate  value  of 
(2ci)  for  the.  truncated  problem,  and  it  was  used  in  Li[*]> 

2)  It  follows  from  (30c)  that  an  exact  solution  of 

LjCw]  •  f2  (A. 6) 

Is  W  *  2f/3.  Let  [Q2n3  denote  the  coefficients  obtained  by  solving  (A. 6) 
approximately,  and  measure  the  error  in  this  approximate  solution  by 

e-  •  max  l*2n  “  2W3I  1  lf0l  (A. 7) 

0  <  n  <  N 

For  N  «  7,  vie  find  e2  <  10-2  for  tr/A  <  X  <  */2,  and  e2  <  6  *  10-2  for 
w/10  <  X  <  3ir/4. 

3)  Truncation  of  the  series  at  N  requires  that  the  first  few  terms  dominate 
the  series.  Let  t?2n  denote  the  coefficients  in  the  solution  of  (A. 2),  and 
define 


e3  »  max  {|v12|,  |v14|)/|vQ|  (A.8) 

For  N  »  7,  e3  <  2  *  10-2  if  v/4  <  X  <  «/2,  and  e3  <  7  *  10-2  if 
n/10  <  x  <  3»/4. 
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